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ABSTRACT 


The need to minimise the pioduction, inventory and tiansportation costs has been long 
recognised by operations reseaicheis I raditionally, the stages in logistical functions 
production at production site, transposition from production site to retailer site and 
consumption at i etailei site have been optimised separately In the present research an 
integrated appioach towards the above functions has been made 

We considei production distubution system with a central production centre and multiple 
retaileis Hithei one single pioduct or multiple pioduct-one for each retailer is being supplied 
through capacitated vehicles The demand is assumed to be continuous and deteiministic The 
objective is to deteimine the optimal leplemshment policy so as to maximise the overall profit 
by exploiting the cost Unde-offs between individual logistical components 

Two leplemshment stialegies single replenishment (i e one replenishment for a retailer in a 
single production cycle) and multiple replenishment have been analysed The focus of the 
present work is on analytical approach to solving production distribution problems The 
solution methodology follows the same line in both the cases It is first decomposed into 
retaileiwise subpioblem and then solved independently Finally these subproblems are 
integrated to foim the simplified objective function The objective function is analysed for 3 k 
ranges of the domain of n (numbet of pioduction cycle) where k is the number of retailers In 
each langc the local optimum can be obtained in polynomial time Then the best solution is 
obtained The computational complexity of both the algonthms are O(^) subject to the 
condition that squat e loot of a number is obtained in constant time Results of a computational 
study on landomly selected piobletns are also presented 



Chapter 1 


INTRODUCTION, LITERATURE SURVEY AND 
PROBLEM SETTING 

1.1 Introduction 

With the trend towards greater syneigy between suppliers and industrial customers, most 
manufactuimg entei puses are oiganised as networks of manufacturing and distribution sites 
that pioccss the raw matenal and distribute the finished goods to customers Management of 
such netwoiks (also iefened to as "supply chains") has recently emerged as a major topic in 
Opeialions Research (Lee and Billington, 1993) 

Optimal management of such networks requnes caieful analysis of the combined logistics 
pioblem The integrated model is used to study cost ti ade-offs between individual logistical 
components vi 7 production, tianspoitation and tetail entity, as well as interrelations between 
vatious logistical functions 1 his study can be utilised 

i) to cieate an accurate representation of trade-offs encountered and 
n) to solve the combined problem 

This type of combined logistic problem will be veiy useful for organizations, that own both the 
pioduction and retail outlet facilities and thus are lesponsible for the replenishments and 
inventory management at both the sites The present reseal ch concerns the management of one 
segment of a supply chain by such an integiated approach decisions for both inventory and 
tianspoitation functions aic made simultaneously It will help in detci mining the pioduction 
pattern and the shipping stiategy for each lelailei foi maximum piofit 



1.2 Literature survey 


Although theie is substantial hteiature on inventory and transportation management 
respectively, much less is available on the combined problem Initial study in this area began in 
mid-1980s 

hcdeigiuen and Zipkin (1984) mtegiated the allocation and routing problems within the same 
mode! by ticalmg a single period, single item pioblem with random demands at the retailers 
No ordeung cost is considered r lheir method is to first decompose the problem into two 
poitions mventoiy allocation and vehicle routing The solution then is obtained by 
constructive and/or intci change heunstics 

Bums el al (1985) developed a single item model with the objective of minimising 
transpoitation and inventoiy cost pei unit tune Production related costs are not included The 
tianspoitation cost consists of a fixed chaige pei unit movement It is assumed that production 
is not synchronised with delivery and the cost of inventory accumulating before delivery is 
estimated accoidingly 4 he optimal delivery quantity is obtained by an economic order 
quantity type of analysis 

The cost tiade-ofTs between vauous logistical components on fi eight netwoiks have been 
analyzed by Daganzo (1985) Networks involving direct shipping, shipping via consolidation 

i 

teiminal and a combination of teimmal and direct shipping are considered by him He 
developed a simple optimizing method that simultaneously determines optimal routes and , 
shipment sizes for netwoiks He also identified conditions that indicate when networks 

r 

involving dnect shipments between many ongms and destinations can be considered on a link I 

i 

by link basis 

t 

i 

[ 

i 



Benjamin (1989) investigated a single item veision ol the economic lot and dehveiy scheduling 
pioblem in which multiple dehvenes may be made out of one batch He excluded from his 
objective function the cost of mventoiy accumulation befoie delivery and the additional 
inventoiies that would be lequired to avoid shoitages when pioduction interval (or batch) is 
not an integer multiple of the delivery intei val (batch) These two implicit assumptions lead 
Benjamin to the conclusion that the pioblem may be solved optimally by independent EOQ 
type foimulas for the pioduction and delivery batches Unfortunately, the independent 
solutions generally will not have the mtegei multiple property that is assumed in the 
fot mutation 

T W Chem (1991) dealt with profit onented oigamsations that own both production and retail 
entities and face demand fluctuations He consideied the case of stochastic demands and 
follows a probability distribution Ills teseaich extends the classic Newsboy model by 
combining the pioduction and logistical decisions He treated the production rate and shipment 
size as decision vanables and calculated expected profits Monte Carlo experiments and 
sensitivity analysis aie also perfoimed to validate the model and view further insights into the 
effects of changes in parametei values on expected profit function 

Anily and Fedeigiucn (1990) considei a single item, deteiministic model where the objective is 
to minimise avciage tt anspoi tation and inventory costs over an infinite horizon They use a 
policy where a retailei can be assigned to different loutes simultaneously, each route satisfying 
possibly a fi action of the total demand A paitilion method is used to deuve lower bounds on 
the lota! cost It is then shown that these bounds become asymptotically tight when the 
numbei of demand points tend to infinity Anily and Fedeigiuen (1993) extend this woik to a 
two echelon dislt ibution pioblem, wheie the cential waiehouse keeps system stocks instead of 


being a meie twmshipmeut point Gallego ahd SimUn-Levi (1990) piesent a lowei bound on 
the long iun aveiage cost ovei all inventory routing strategies for such systems They stress 
the benefits of dnect shipping, wheieas Hall (1992) atgues that multiple-stop i outing strategies 
can still pi ovule laige savings 

For the case ot deteinumstic demands foi several products at several retailers, Vishwanathan 
and Mathur (1994) consideied both inventory decisions and vehicle routes from a warehouse 
They obtain a stationary nested joint replenishment policy for the situation when replenishment 
mtei vals are powei of two multiples of a planning penod 

V Anand (1997) has consideied the combined logistics pioblem for a production distribution 
system with single production centre, single letailer outlet and a unique product under 
detemiimstic conditions He has formulated the problem and has derived closed form solution 
under uniform as well as non-uniform shipment policy 

1.3 Problem Setting 

In the present leseatch, we consrdei the profit oriented organizations, which own both the 
pioduction and telail entities The goal of the present research is to determine the optimal 
shipping strategies which maximises the overall profit by analysing the trade-offs that exist 
between production, inventoiy and tiansportation costs To obtain shipment policy, we 
foimulate mtegiated model that consists of all the logistical components 

'Hie pioblem consists of pioducing a unique product (or multiple product, one for each 
ictailci) at production site and shipping it diiectly to the vniious retailer using capacitated 


vehicles I he demand at the outlet is assumed to he umlmuos mid deleimmistie in case of 


multiple pioduct, it is assumed that there is one identified letailer for that product 

At the ptoduction site, the items aie produced at a constant production rate The shipment to a 
retailei site is of equal volume Howevei the volume and number of shipment will differ from 
retailei to letailer The costs incuired at the production site are cost of production set up and 
the cost of inventoiy holding 

The capacity of transportation vehicle is limited There exists a fixed transportation cost for a 
shipment and a vaiiable cost propoitional to the number of units transported These charges 
aie retailer dependent An inventory holding cost for the duration of transportation is also 
lncuited 

A tegulai revenue is obtained per unit of demand satisfied at the retailer outlet The demand 
not satisfied is consideied lost and theie exists a pei unit shortage cost The inventory holding 
cost at the letail outlet is also taken in account 


1.4 Organisation of thesis 

In chaptei 2, the case when only one leplemshment is made to retailer site in a single 
production tun is consideied In chaptei 3, the pioblem, m which more than one replenishment 
can be made to a letailer site in a production run is analysed In both the chapters the 
mtegialcd pioblem is lot initiated, the constraints aie identified and then a methodology is 
developed to solve the problem In chapter 4, we btiefly discuss the conclusions and 
suggestions foi futuie vvoik 


Chapter 2 


PRODUCTION DISTRIBUTION SYSTEM WITH SINGLE 
REPLENISHMENT 

In this chapter, we considei production distubution system with one production centre and 
multiple letaileis r !he shipment fiom the piocluction centre to the retail centre is made 
through vehicles (ships) with limited carrying capacity We further assume that only one 
shipment is made in each production cycle We consider the problem of finding the 
optimal shipment policy i e number of shipments per unit period, and the volume of 
shipment for each letailei, so as to maximise the oveiail profit aflei taking inventory, 
ptoduction and tianspoitation costs into considerations 

In section 2 1 we determine the general relationship between production pattern and 
shipment policy 

2.1 Production Pattern and Shipment Policy 

The system consists of a ptoduction site, a numbei of retaileis and one capacitated vehicle 
foi each retailer 


Retailers 

► 

*• 

Figuie2.l A Simple Production Distribution System 
2.1.1 Assumptions 

We consider this case with the following assumptions 
1 Demand at a letailei is continues and uniform 


Production 

Centie 




2 Production late is greater than the cumulative demand rate of all the retailers 

3 Demand not satisfied is lost 

4 Net revenue is always positive 

5 Only one shipment is made to a retailer m one set-up 

6 One vehicle is available for shipment to each letailer, le no scheduling of ships is 
involved [In section 2 8 we examine the case wheie the number of vehicles is not a 
constiamt] 

7 System is pei petual 

We have considered the following cost and revenue for the purpose of optimisation as 
they aie likely to be afiected by the shipment decisions 

A) Pioduction site 

i) Production set up cost 
n) Inventoiy holding cost 

B) Transpoitation 

i) Cost of tiansportation 

u) Cost of transpoitation inventory 

C) Retail outlet 

l) Revenue leceived 
n) Cost of holding inventoiy 
ni) Cost of lost demand 

2 12 Notation 

I he following notation is used 

Pi obfent specification data 

Xj Demand rate per unit tune period of retailer j 

P Production capacity of the production centre per unit time period 

IV, Maximum shipment capacity of the vehicle for retailer j 

R pei unit net revenue (gi oss revenue- cost excluding inventoiy and 

ti anspoi tation cost ) 

T, 1 1 anspoi tation time foi ictailci / 



C Fixed cost pei production set up 

bj,jj Fixed and per unit variable transportation cost of shipping to 

retailei j 

s per unit shortage cost 

hp, lit pet unit holding cost pet unit time period at production centre and 

during transportation respectively 

hi j per unit holding cpst per unit time period at retail outlet j 

k total number of retailers 

Decision variah les 


V) Shipment size foi retailei / 

yopt 

Optimal shipment size for retailer j 
n number of set-ups in the given period 

Hops Optimal value of n for the given lange 

6j binaiy vanable , 0 if replenishment is not made to retailer j 

1 if replenishment is made to retailer j 

2.1.3 Observations 

When the system is peipetual, the following observation is made 

Lemma (2.1) In a perpetual system, all cycles will be of equal length and total profit 
Pi oof Consider one such system in which cycles are not of equal profit for a 

finite peiiod of time 7T Let it consist of // distinct cycles each with profit 
obtained per unit tune is C\, C 2 , . C» Let the length of the cycles be Tj, 

r 2 , . 7’ 

n 

Total profit for the finite length of time L is ^C t T t -(21) 

/=i 

n 

Total profit per unit time period is (]T C] 7 j)/L (2 2) 

M 

Iheie exists a maximum piofit incut ring cycle with profit C max as max [Cj] 
Then 'Zap; < < cW £ 7 I ^ Crtiax L .(2 3) 

l-l I~1 1=1 



=>(£(; ?; )/ls (2 4) 

1=1 

To maximise the total profit the system will have identical cycles of profit 
( T 

Hence the peipetual systems will have identical cycles 


Let us considei the system for a unit time period in which demands at the retailer outlet j 
is Xj and the pioduction capacity is P 'I he shipment of size Vj is carried from the 
production site to the retail outlet j in each cycle, and there are n such cycles per unit 
penod The total amount shipped to a retailer is less than or equal to the demand of the 
tetailer in the unit period 

With the above framework of the pioblein, the following observations are made 
Lemma (2.2) When single shipment is made in a pioduction cycle to each retailer then it 
is always advantageous to produce at the maximum rate for each retailer 
Proof Let the production be done at the rate qj for the retailer j subject to the 

const! aint that o) is lower than the pioduction capacity P In that case we 
will have inventory holding cost at the production site 


= hp 



2 a , 


'1 he conesponding cost at full pioduction rate P is hp 


V 1 
2 P 


As a, <P hp — > hp— 

1 2 a, IP 

Hence the minimum inventory will be when o) ~ P 



Figuie 2.2 : Inventoiy Status Diagiam at Production Centre for both Policies 



Hence foi the optimal shipment policy, the pioduct will be nianufactuicd at the maximum 
possible pioduction late 


Lemma (2.3) In the optimal shipment policy, shipment to a retailer will be made as soon 
as the pioduct has been pi educed (at the maximum rate of production) 

Pi oof Let Vj be the optimal volume to the /' letader Let the shipment be made at 

an mtei val i subject to the condition that V/P < i The inventory status 
diagram lot this case is shown m fig 2 3 



t > 

Figuie 2*3 : Inventoiy Status Diagram for Independent Schedules 


With the same letailei inventory status, same volumes and same number of 
cycles the coi responding inventoiy status diagram at the production site for 
the condition when shipment is made at made time /* V/P is shown in fig 
2,4. 




Figuie 2 4 : Inventory Status Diagram for Synchronised Schedules 

The area of the inventory status diagram gives the inventory holding cost It 
is evident fiom the figures that the inventory holding cost is lesser in the 
second case as the area bounded between V and i is lesser than in first case, 
while all the other costs aie same 

Hence the shipping policy that ships the required quantity as soon as it is produced is 
always better than delayed shipment 

In this section we develop the total piofit function for such a policy 

2.2 Problem Formulation 

In this section we foimulate the net profit function per unit time period The costs and 
revenues at pioduction site, letail outlet and during transportation are calculated 
sepai ately 

A) Production t elated cost 

The following figuie lepresents the model at the production site The relevant 
costs aie cost of pioduction set up and inventory costs 



V 



Figuic 2.5 : Inventory Status Diagram at Production Site 


1) Production set up cost There are //, production set ups in a unit time As C is the cost 
per pioduction set up, pioduction set up cost per unit time is nC (2 5) 


2) lnventoiy holding cost T he total inventory holding cost in one period can be 
computed by measuring the area of figure 2 5 

1 he total inventory cost per unit period is 


~ — !ML V |/ 2 
2 P } 


(2 6 ) 


B) Tiansyoitaiion related cost 


1) Cost of tianspoilation 

A fixed freight chaige P s and a variable cost f, pei unit transported is 
incun ed foi the shipment to retailer j Thus the total transportation cost is 

h tf> v j ] < 27 ) 

r \ J=l 

wheie Oj is a binary vanable, 

= 1 when V } > 0 , 

=0 when V } - 0 

2) Tiansportation inventory holding cost 

If the transpoitation time is 1) for the /' retailer then the total inventory holding 
cost duung ti anspoi tation ovei unit time penod is 

nhtj^VjT, . (2 8) 

j i 


C) Ret ml a outlet 


1) Revenue 

k 

The total amount of demand satisfied in one period is nYV The total revenue 

k 

generated is n R ^ V I (2 9) 

2) Inventory cost r I he following llguie represents the stock position at any time 
at the retailei outlet 



Figuie 2.6 : Inventory Status Diagram at Retailer Site 


If (he per unit holding cost pei unit time period at the retailer site j is hr } then the total 
inventoty holding cost pel unit time penod is 


- Y 
2 




2\ 


X 


(2 10 ) 


, J 


3) Shoitage cost 

The number of lost sales in one penod at a retailer j is (X r nVj) If the per unit 
shoitage cost is ,s then the shoitage cost pei unit time period is 

-sj^(X,-nV 1 ) . (2 11) 

; i 


The profit function is obtained by subtracting all costs from the revenue obtained The 
ptofit is a function of // and K/s Let us represent this function by Z Then 

» r ’> • " fi t",- -riftr; - »('■; t», ■ W -»-'±v,r, 


r i 








I Juf hi.v ^ 1 




J J 

1 I V X, J 


J I 


(2 12 ) 







2 \ 


I V 


X. 


;=l jn 1 ;=1 



(2 13) 


We will maximise Z to obtam the optimal policy The feasible policy has to further satisfy 
the following const! aints, 

1) Supply is less than oi equal to demand 

nVjZXj (2 14) 

2) I he volume shipped is less than or equal to the vehicle capacity 

V } < Wj - (2 15) 


3) Only one vehicle is allotted per retailei so that production cycle must be so 
spaced that the vehicle can return 1 e 


n ^ mm[ — ] 

e>r, 

(2 16) 

4) non negativity constraint 


Vj> 0,n^0 

(2 17) 

2.3 Solution Methodology 



Optimization of Z includes decision variables n and V } We solve the problem by first 
getting optimal Vj Jpl for each retailer as a function of n and then optimising it over the 
decision vanabies n to obtain optimal values i e n opt and V° p 

For a selected value of n tD, the pioblem reduces to maximisation of Zj(V) We refer to 
this pioblem as A 


Problem A 


;=) L1 1 


k y 2 

j 

P 


Z 




1 V 


X, 


-hi±vn±FA-±tV,-C] 

j* 1 ;=! j= 1 


subject to the const! alula as given by Eqn 2 14, 2 15 and 2 17 


.. .(2 18) 



It can be obseived that objective function can be written as Z(V) = y^Z(K t ) Further 

there aie no constiamt binding the different decision variables together so the problem can 
be decomposed in solving independent subproblems each of the given type 


Subproblem SP 

The reduced f' subpioblem is 

Z 2 = (R i s-fj) Vj 


!hj_ 

2 P 1 " 2X, 


= (R<s-f r htT J )V J -(^ + 


V/-IU V/lfF&j 

hi , 

2X>*- f 4 


Jip hi 

Let a," + — — and $= R t s~f r htl) 

2P 2 Xj 

r I hen, substituting a } and p } we obtain 


(2 19) 

(2 20 ) 


Zj- Wr^Vi-hB, 


(2 21 ) 


2 3.1 Observations 

Case 1 If p } ^ 0, then V } =-0, otherwise the losses met ease on increasing replenishment 
size 


Case! If/? >0 then £ _ F K '>0 (2 22) 


As we are now concerned with only one ictailei we can drop the subscripts 
From the equation 2 22 it is cleai that only those values of V have to considered for which 
7.2 >0, because if lot a V ' 0, Z 2 " 0, then it implies that total cost (inventory holding cost 
and tianspoitation cost) is highei than net tevenue at that volume of shipment, and it will 
not be viable to make shipment at that values as V~Q is a better solution 

Next we will calculate the range for V which Z 2 > 0 by finding the root of equation 


z 2 -o 


Foi Z 2 to positive, V should take values between 


p ± yjp 2 - AaF 


2 a 


(2.23) 


If jf ~ 4 ah' < 0 then the roots aie imaginary which implies that then Z 2 can not take 
positive value , so leplenishment will not be made in that condition 


2 . 3.2 Solution method forSP 


One can look at the nature of the function which can be either of the two types as given 
below 



Figure 2.7 : Variation of Z with V. 


From figure 2 7 it is clear that V Qpt - 0 for case 1 Now we will consider case II 
From unconstrained optimisation of Eqif2 21 we get 

= p-2aV . (2.24) 

dV 

< 0) so the function is concave and the K* corresponds 

to the maxima of the function 

B 

So, the unconstrained optimal V is V* (2.25) 

2 a 

Thus the value of permissible V varies with n as given below in the figure 


d 2 Z a 
dV 2 


= -2 a ( 



Figui e 2.8 : Representation of Permissible Range of V for any Retailer 



I he value of Lis bounded by W, the vehicle capacity on the uppei side and by Vr ouer on the 
lowei level 

A f X \ 

Let V =111111 l ' — ,W then h om the above analysis it is evident that 
\ n J 

V/ ,pt = V if V > V/oner and n > 0 

— 0 V < V/ aite r or n—0 (2 26) 

1 he figure representing Vj" pl for different values n is as given below 


♦ 

v 



a). When V* > W 



n 


b) When V* < W 

Figure 2.9 : Dmgiam Showing V t) | )( at any value of n 


So, we have at most two bieakpoints for each retailei with respect to n The first shows 
the point upto which the replenishment is of constant level (independent of n) min (K* 
ff) 1 he legion between the first breakpoint and the second breakpoint gives those values 



of n foi which full teplenishmcnl ( Y'/t) is to be made Beyond the second bieakpolnt the 
size of optimal icplcnisimicnt is zcio 


( 


I fie first bieakpoint, B\, is achieved at min 


1 




(2 27) 




I he second bieakpoint, B 2 , is absent when — < - 

7’ min^V^j 


otherwise it is attained at 


mm 


* 


w’r ) 

Heie l r i imt ,r coriesponds to the lowei of the value obtained in Eqn 2 23 
I hus at this stage we can compute Vj° p ' of any letailer j for any value of n 


(2 28) 


2.3 3 Solution method for Problem A 

The retaileis foi which optimal solution for all values of n is no shipment (i e V J opt ~0') need 
not be considered any luither Let G lepiesent the set of index of all such retailers 

Obiective function : 

Foi any given //, the objective function can be written as follows, 

z4'd'-z z ^> < 229 > 

;cG, jri tj 

whei e 

Gi^ { all / such that Vj=. min(K* fF) }and 
G 2 = {all / such that Vj^X /11 ) 

As at most theie ate two ranges for each retailei, the whole range of n (0 to 00 ) can be 
petitioned into at most 2/r+l distinct range For each, the objective function can be 
calculated by identifying the set G| and G 2 I he objective function foi a specific range will 
be one of the following type 

A) Type I G 2 ={0} 

When all ictaileis aie being supplied at a level which is constant (independent of n )then 
the value of the function Z\ is as follows, 

k 

Z ~ /j[ constant] - (2 30) 

j 1 «. 


Z is a linear function of n and hence within this range, n would take the upper limit of the 
uuige as the optimal value 'I hus foi lias mnge the local optimal value is the uppci limit of 
the tango 


13) Type II : G i“{®} i a all i etailers aie being supplied at full demand level (those having 
no replenishment have been all eddy neglected) 


X 

Substituting Vj = — in Eqn 2 13 we obtain 
n 

z 4 („) ^±X,-^±X)-±±h, 1 X l -±$X,-«±F,4a±X l T l -C 

)~\ j~ I 1,1 j=\ j*\ ;=1 


where j <£ G 


(2 31) 


Let wt* - tfX -totx.h ■ e^±X] t ±±ln t x, and ft«£i5+C 

j I z ;=l ;=l 


) i 


}-- ) 


r ' 


where j & G 


then the above equation can be expiessed as 

z<(n) = 

whete 0? and O 3 aie non-negative numbers 
By differentiating Z 4 with respect to >?, we obtain n* - 


\6i 


(2 32) 


.(2 33) 


and as the second derivative ( — -) is always negative the given function is concave 

n 

Now we check the range for which full replenishment is being made 

If n* > upper limit of the range then n op , = upper limit , 

if //* < lower limit of the range then n 0])l = lower limit else n opt - n* 


C) Type III : Gj and G 2 ^<i> / e Some retailers are being supplied at constant level 
(independent of «)and some aie at full replenishment level 
In this case the Eqn 2 14 can be written as 

k 

Z - nlZi(Vj.) ]- ij^x, where /* eO,, G 2 (2 34) 

} 5 

The contribution due to set Gi is a lmeat function with lespect to n while the contribution 
due to set G 2 is of the nature as discussed 111 the above case (i e Z^n) - 6i~(9 2 /n) -w#j) 
Let the contubution from set Gi be n9r$s then we have 
Z- n(QrGs)-(6 2 m) ' 9, -0 S 



as in the above case 

If $ 4 ~ 0 s > 0 then the function always mci eases with n , so in this case upper end of the 
tange is chosen as then 0}Jl 

If OrOh < 0 then on the same pallem as in type II, we obtain 


//* = 


Oi 


03 ~ 04 


(2 35) 


Now we check the tange foi which full replenishment is being made 

If //* > upper limit of the range then n op{ = upper limit , 

if//* < lower limit of the range then n opt - lowei limit else n opt - n* 


2.4 Algorithm 

Fust we independently solve the letaileiwise subpioblem and then use this result to solve 
the total problem 1 he detailed algorithm is given below 


Subproblem SP : 

Initialise / ~ 1, G={<D} 


Jjn 111 

a) Calculate a } and ft } wlic/e a } ~ - - p — and ft } ~ R i s~f r nil) 

2P 2Xj 

b) If ft } < 0 then assign j to set G and go to step g 


c) Calculate v s(ai ft] 


_ p, 


2 a 


and v lowei j]] 


_ PjzJKzXA 


2 a 


if ft , 7 4 (Xj h) then assign / to set G and go to step g 

d) If v slat ft] Wj then v slai ft] W } , 

e) If r s uv ft] v lowei ft] then assign / in set G and go to step g 

f) Calculate the bieakpoints as follows 
, X. 

Rift] - mm 


T } 

if Bl[j] then B2{j] = B 1 [j] elseB2[j] = min 

^ i 


X. 


Tft v lowei [i] 



g) Repeat the procedure for the next retailer 


11 Problem P 


a) Combine at rays For / =1 to k and / 0G combine B 1 [j] and B2[j] and 0 and then 
arrange them m an ascending manner 

b) Make pairs out of the above arianged list from each two consecutive numbers 
to const! uct ranges 

c) For each range, check its type 

Calculate Gi and G2 

wheie Gl= { ally such that Vj=min(y*,lV) } and 
G2= {all / such that V f ~X/n } 

If G|={0} then range type is Type I, 

If <32={<t>} then range type ts Type II else range type is Type III 

d) Calculate n opt for that range according to the range type 

1) 1 ype I 7 t op , = langejippei 


11) Type II n opt = 


- P ~Y x]+~Yh, x, 

o T) j O ' 7 7 


2 P% 


Z 7=) 


where j 0G 


7= 1 


if 11 opt > wnge uppet thep n op , = icmgejippet , 
if }i 0 pt <• / angejowei then n opt = 7 cmgejower 

iii)Type III Calculate the conti ibution to objective function due to set Gi 

jt 

I ?))0 4 -&s is the contnbution, then if 0 i> , where j 0G then 

/=! 


7 i opt = langejippei 


else 


Calculate Or*~Y, x ] + ^'L h, i Xj 

II j-i 


9 x — 
z 7=1 


Hope 


e. 


whereyeG2 then 
where j 0 G 


7^1 


if 11 op t > iange_uppei then ir op , = range_upper, 
if 7/ w/J ( < rangejower then n op , = rangejower 

e) At each of these n 0 pt the cot responding F/s aie calculated 



if n ol n < B 1 [)] then V? A - min(K*, W } ) 

X 

if n ol >t lies between Bl[j] and B2|jJ then VJ >vt = — ~ 
if IV > B2 (j] then F/ /,l =0 

f) At each n opl calculate the objective function value and select 
which objective lunction value is maximum 

2.5 Flow Chart 



Figmc 2.10 : Flow Chait for the Algorithm 









2.6 Complexity Analysis 


Let k indicate the number of letailers In the first part, SP, each of the five step requires a 
constant number of operations (addition, subtraction, comparison, division etc ) so the 
total number of operations for one retailer is a constant, say L This procedure is repeated 
k numbei of times thus the numbei of opeiations requiied is Lk, which also assumes that 
square root opeiation takes a constant amount of time 1 e operation tune is independent of 
problem size 

Foi the pioblem P we have to first sort the 2k ranges which requires (2k log 2k) steps 
Now for each range we make k evaluations to find the optimal decision variables This 
step requires 21c steps 

Thus the pioblem lequires Lk+2k log 2k \ 2k 2 steps Hence the complexity of the above 
algorithm is 0(k 2 ) i e it is a polynomial function of number of retailers 

2.7 Example 

1 o lllusti ate the above analysis we consider a numei ical example with the following data 
R= $615 per unit, c= $600 per unit, P= 5000 units per year, s= $10 per unit, hp= $80 per 
unit per year, hr = $50 per unit per year (for all retailers), ht = $15 per unit per year, k= 4, 
F,= f 2 = F 4 = $1000 per transportation , f 2 = f 3 = £»= $1 per unit 

Xi= 600 units per year, X 2 = 900 units per year, X 3 = 1200 units per year, X 4 = 1000 units 
per year, Tj=0 05 year, T 2 =0 08 year, T 3 =0 125 year, T 4 =0 10 year Wi= W 2 = W 3 = 
W 4 = : 450 units 

Solution \ As explained in the above section we can get a table like the given below 


Retailer no — > 

1 

2 

3 

4 

v* 

234 

318 6 

383 7 

340 91 

V ]owtr 

47 91 

47 38 

48,22 

47 79 

min (V*,W) 

234 

318 6 

383 7 

340 19 

Eh 

2 56 

2 82 

3 13 

2 93 

b 2 

12 5 

125 

8 

10 

Table 2.1: Su 

jpioblem Calculation 






















On ananging (he combined set of bieakpomts in an ascending senes and then making pair 
out of the adjacent numbeis we get seven ranges The subsequent analysis is presented in 
the table below 


Range 

Type 

n Dp ( 

Z 

0-2 56 

I 

2 55 

-10291 65 

2 56 - 2 82 

III G.={2,3,4} G 2 = t { 1 ) 

2 81 

-7742 

2 82 - 2 93 

IllG l ={3,4},Gr={l,2} 

2 92 

-6665 

2 93 - 3 13 

III G I ={3),G 2 ={ 1,2,4} 

.. 

3 13 

-4863 

3 13 -80 

11 G 2 “{ 1,2,3,4} 

5 44 

2447 

80-100 

II GH 1,2,4} 

8 0 

-14012 

10-125 

11 Ga— { 1*2} 

10 0 

-17994 


Table 2.2 : Pioblem Calculation 

From the above table it is evident that the optimal value of n is 5 44 and the replenishment 
sizes for the retailets are HO 27, 165 41, 220 54 and 183 79 units respectively 

2.8 Problem with Unlimited Number of Vehicles 


7 he same problem that has been considered in the above sections can be solved for the 
case when the number of vehicles is not restricted This implies that the constraint that 

n <> min[ 1 need not be satisfied 

87 ; 

Due to this relaxation, the bieakpoint calculation is altered slightly 

( \r \ 

The first breakpoint, B|, is achieved at 


X 


7 he second bieakpoint, B 2 , is attained at min 


i,min(K*,^) 
X 




(2 36) 
(2 37) 


Ilcie I'/u,,,, couesponds to the lowci of the value obtained in Eqn 2,23 
I he lest of the solution methodology is exactly the same 




Chapter 3 


PRODUCTION DISTRIBUTION SYSTEM WITH MULTIPLE 
REPLENISHMENTS 

In this chapter, we considei the pioduction distribution system where multiple shipments to a 
retailers in a single production cycle are permitted We consider two kinds of production 
policies In the first situation, only one production setup for each retailer, in each cycle is 
permitted and the production takes place at the maximum production rate P Such a situation 
may occur when the production is in identified batches for each retailer (1 e order processing) 
or multiple pioducts are being produced, one for each retailer In the second situation we 
consider that simultaneous production take place for each retailer and shipments are made to 
each of the retailers throughout the cycle Such a case may refer to a situation where a 
common product is being produced to be shipped to all retailers or to the case when all the 
products are being produced simultaneously at the rate consistent with their total 
replenishment m the cycle 

As in the previous chapter, we first determine the general behavior of production and demand 
patterns for the above situations and then foimulate the profit function by considering all 
relevant costs and tevenues Next, we optimise the profit function subject to certain constraint 
to obtain the optimal shipping policy 

Assumptions 

'I he same set of assumptions as in section 2 1 have been considered except the assumption of 
only one shipment per pioduction setup lias been relaxed 

Notation 

In addition to the notation described in section 2 1 we will use the following notation 
m, numbet of replenishment in a single setup for retailer j 

jo if hp it hi, 

1 [1 tf lip < ll) J 



The nature of the shipment policy will depend upon the relative cost of maintaining inventory 
at the production site {hp) and at the retailer site (h/) We consider the following two cases 

Case l If the cost of maintaining inventory at production site is larger than the retailer site 
(1 e hp > hij ) then optimal numbei of shipments of fixed size will be made to the 
letailer centre as soon as possible so as to reduce the inventory at production site at 
the cost of increasing it at the retailer site 

Case II If the cost of maintaining inventory at retailer site is greater than that at the 

production site (i e hp hfj ) then the shipments will be delayed as far as possible so 
as to reduce the mventoiy cost at the retailer site at the expense of increasing inventory 
holding cost at pioduction site 

Next we consider the production policy m which only one production setup for each retailer is 
pei nutted in a single cycle 

3.1. Single Production Setup Multiple Replenishment for each 
Retailer (SPS) 

In this case we consider the situation when total order of a retailer is produced at the 
maximum rate of production only once in a cycle The product is shipped in shipments of 
equal size that aie spaced at equal intervals throughout the cycle We shall formulate the net 
profit function foi this case I he costs and revenues at pioduction site, retail outlet and during 
transportation are calculated separately 

3.1.1 Problem Formulation 

In this section we fbi mutate the net profit function per unit time period The costs and 
revenues at production site, retail outlet and dui mg transportation are calculated separately 

A) Production related cost 

1 he following figures teptesenl the models at the production site foi case I and case II 
j especti vely 1 he i elevant costs ai e cost of production setup and inventory costs 



Case 1 



Figuie 3.1: Inventory Status Diagram at Pioduction site 

1) Pioduction set up cost There are ft production setups in a unit time As the cost of 

set up is C, the total production set up cost per unit time is nC (3 1) 

2) Inventory holding cost The total inventory holding cost in one period can be 
computed by calculating the aiea of the above figuie If the inventory holding cost is hp then 
the total inventory cost per unit period for case I is, 

. 1 ‘ m,Vf 


2 ^ P 

l 

11 hp ^ r/2 

= — 

2 P Yl } J 


(3 2a) 


and for case II is 


///?/;* K 2 / x ( ill X \ 

Z m > y- ("'/ + 0 - 1 ■ -y-j ISce Appendix-A 1 ] 
j L 


(3 2b) 


i 


! 


li) Ti ansportation related cost 

1) Cost of transpoi Cation 

A fixed freight charge F } and a variable cost f } per unit transported is incurred 
depending upon the retaiiei site I hus the total ti ansportation cost is 

h 1 ' ( 3 3 ) 
j=i j=] 

2) Ti ansportation inventoiy holding cost 

If the transportation time is 1) for the j lh retailer and the transportation 
inventory holding cost is hi then the total inventory holding cost during transportation over the 
unit time period is 

k 

nht Y.nijVJj . (3 4) 

j = i 

C) Revenue and cost at the retailer site 

1) Revenue 

k 

The total amount of demand satisfied in one period is nijV 1 The total 

j=\ 

ievenue generated is 

nR±n,V } (3 5) 

2) Inventoiy cost 

The following figures represents the stock position at any time at the retailer 
outlet for case i and case II respectively 





Figure 3.2 .• Inventory Status Diagram at Retailer Site 


If the per unit holding cost at the retailer's site is hij, the total inventory holding cost per unit 
time peiiod is 

for case I 

« ^ ) [See Appendix-A 2] 


and for case II 


1 * h t ttbV 2 
n — > ~ 

L )=\ 


(3 6) 


3) Shoitagc cost The number of lost sales in one period at a retailer ; is (X r ti WjVj) 
If the per unit shortage cost is s then the shortage cost pel unit time period is 

s jr(x,-mV,) ' (3 7) 

;= « 


The piofit function is obtained by subti acting all costs from the revenue obtained The profit is 
a function of ti, m } and V/s Let us lepieseut this function by Z t Then we describe this 
pioblem, denoting it as P, in the next section 


3.2 Problem P 


Z,(n, in, Vj) = nR J^niiV, - nC ^ Z( 1_ fO"'// 


2 P 
f /;/j /??, - A ;i/?p 


Z ;=1 




* r/^ 

2 


h +l )-Kr- 


| * (h hi IlhV 2 k k k k 

— --"[Z^ ^ " lu ^L n ') V ‘ T i - s^(x - imijV,) 

1 ) \ X ) j-\ 1 j=l 


(3 8) 


We will maximise Z; to obtain the optimal policy The feasible policy has to further satisfy the 
following constraints 

1) Supply is less than or equal to demand 

n ntj V] <X } (3 9) 

2) The volume shipped is less than or equal to the vehicle capacity 

(3 10) 

3) Only one vehicle is allotted per retailer so the production cycle must be so spaced 
that the vehicle can return to take the new replenishment 


mtj ^ 


1\ 


4) non negativity constraint 
ii, nij Vj> 0 

The function Zj(n, m } Vj ) can be rewritten as follows 

Zi(„, m, !',)=» {R imV, - C 

}=) )=\ 


1 ^ 

2 






; i 


r/ 2 r 

j 

X. 


{’”i + 0 


ill, X t 

7~p~ ) 


(3 11) 


(3 12) 


i * 6 hi m,y 2 it i k k 

{z— 7— -cz^ 1 Z'"^ i - ,u 'L m > v > r ‘ + s Z w ^ >- *Z x , 

" j i y * j 


(3 13) 


f i 


j i 




j* i 




Thus Z;(^^^)=u2Z M (m ; ,Fj 



^ to- #>) -nr)- 


/?/; 


F 


t r 


'///.A', 




1 * ^/ V »,F, 2 * 


2 ft’ ' ' ,r ; i; '■ ’ K P ) } 2% X, 

AmJ £~J ' j u' ' 

j=\ /=! 


(3 14) 

subject to 


1 )MjVj*X/n 

(3 15) 

2)0 si] <W 

(3 16) 

3) 0 s m, < — t- 
»7, 

(3 17) 


Further m this equation as no constraint binds different decision variables together so the 
problem can be decomposed into k independent problem, one for each retailer 


3 2 1 Subproblem SP 

Zn(m„ V) - (R+s-f r to TJ m,V, - (l ■ - <t >, ) (% + 


hp.hifm 2 


2 P 2 \Xi 




hp + hi j 


,{ p - x ,) 

IPX. 


+ 1 




X . 


- F J»b 


V 2Jf, J| 

(* i s-f r hn)) nijVj-(l - f), ) + ^) - Fy», - 


(l -{/>,) m 2 V 2 hr 


P-X. 

k ~2PF } 


/ 


hp + ht 
2X. 






V 


IPX 


/ ; 


(3.18) 


Let a, = (l - <j> , ) + + , ^' y y - j , P/= 7? +s-jJ- to 7) , and 


Tj 


=(>- 0 | 


r / j -,o 


to, + ?>, 


P- 


y IPX } j 


, 2 AT, y 
It may be noted that ctj and yj > 0 


/ Ip 


On substituting a j5 Pj and yj we obtain the following expression 
ZuOiij, VJ = P, m 1 V r a J w ] Vj 1 -m l 1 Vj'‘y l - Fy», 


(3 J9) 



3.3 Solution Methodology 


As the objective function is separable in k subpioblems for a fixed value of n, we will solve the 
problem as follows 

i) we will study the natuie of optimal solution for the subproblem SP for a fixed value 

of/; 

11 ) we will optimise the objective function with respect to the decision variable n 

3.3 1 Observations 
Case A 

If pj ^ 0, then nij and Vj OI>, ~ 0, as for any V } > 0 the profit will be less than zero 
y"i Jt =o :=> m = 0 

Case B 

If pj > 0 then the equation 3 19 reduces to 

Zi 2 (ni J} VJ - 0 for V } = 0 

m Pj tiijVj-ajmjVj 2 ~Wj 2 Vj 2 yj - F } mj for V } > 0 (3 20) 

If 1)^0 then mp ,( = 0 

Fuither if V } > 0, then there will be atleast one shipment per cycle and hence m } ^ 1 

Optimisation for a Fixed Value of n 

As we aie now concerned with only one letailer we can drop the subscripts and will consider 
the function 

j pmV - amV 2 - m 2 V 2 y - wF V > 0,/» £ 1 
w/w. ;-| 0 K = 0,w=:0 

Next we shall examine this function in the domain V> 0, m £ l 

Fust we will calculate all the points which satisfy the fust order necessary condition for any 
point to be maximum Among the points so obtained, we will compute the set of local 
maximum points, if any 'I he fust older necessary condition for a local maxima (or mimina) 


are 


= pm - 2am V - 2 ym 2 V = 0 
dV 


(3 21) 



(3 22) 


Hr 2 * =PV- aV 2 - 2)i)i V 2 - /< =» 0 

As we are consideung only the legion where V s 0 and m > 0 (m > 1) we obtain from eqtvklrs^ 
3 21 the point which may be a local maxima (or mimima ) iff 

i/* _ P m 


2(a + pii) 

Fui ther substituting it in eqn 3 22 we obtain 


(3 23) 


'p z a B 2 a 

a k a 

4F V 4F 


(3 24) 


As we interested only in the points where V ^ 0, m > 0, there is at most one point 


— j?) 

y 2(a + pw) 

which has the potential for a maxima 
Furthei substituting /;/* in eqn 3 23 we obtain 

I p 2 a 

(,„*p*)=(XM — _ l£) 

y V « 

The //?* and V* should additionally satisfy the constraints that V*< W and m*V*<X/n 
So the region where m* and F* can lie is 


(3 25) 



Figure 3,3 : Feasible region of m* and V* 



Now if Z/ 2 (/»* V*) > 0 and it also lies within the intenor tegion, then the point (w* V*) gives 
the maximum, because at m 1 and K=0 the value of Zi 2 < 0 and, there is only one point 
(w * K*) which satisfies the necessaiy condition for a local maximum or minimum 
Hence we conclude that the optimum solution will be one of the following type 

a) V ~ 0, m = 0 if P < 0 or Z 12 (m* V*) < 0 

b) (m*,V*) ifZi 2 (m^V*) > 0 and (m*,V*) lies within the boundary 

c) At the boundary (if it lies outside the boundary ) 


3 3.2 Search at the Boundary 

1) At the boundary m= 1 

From eqn 3 23 we get V*- ~- — — - (3 26) 

2 (a + r) 

but V >vl should also he between W, the vehicle capacity, and the oversupply 

const! amt i e V° pt - min( • - — — r, W, X/n) 

2 {a + r) 


2) At the boundaiy V ~ W \ m > I and mV <X/n 
Substituting V- W in eqn 3 22 we obtain 

'*_ f3W-aW 2 -F 

2 yW 2 ' 


m 


if {fiW-aW 2 -}'') < 0 then /»*=! as m > 1 
if [flW-aW 2 -F) > 0 and if m*W-> X/n then 


so m 


op l _ 




1 

piv 

- aW 2 - b 


2 w 2 y 

X 


nW 



pW~o\V 2 - F <0 

pW-aW 1 -FZO 

PW - aW 2 - F X 
2 yW n 


3) At the boundaiy mV X*n 

We set mV-- X/n in Eqn 3 21 to get 


(3 27) 



7 px aX 2 yx 2 

Z\ 2 = “ ; — - mb 


ii n“m n 


m 


12 


aX 2 
dn wV 


-/! 


<? 2 Z 


12 




(3 28) 


, (3 29) 


dn irm 
Equating the first derivative to zero we get 

X la 


( < 0) Hence the function is concave 


m 


*=. 


nV* 




if V* > W then V* = W and m * = XI{nW) 

So m" 1 " = min ( — M ,A7 QiHO) and K"" = mm J W 
n \ h V V o: 

1 


4) At the boundary m 


tiT. 


On substituting this value of m in eqn 3 23 we obtain 
( \ 


V * - min 


Xl\ 


npr, 


V 


(3 30) 


■"2 (nal] + r )j 

Tims theie are at most six points which can be the candidate for the maximum, four on the 
boundaries, one at origin and the last being the unconstrained optimal 



Figure 3.4 : Candidate Points for Maximum in the Feasible Region 



Al these points, the vauous distinct values which m* and V* can take is given below The 
numbei in bracket indicates the coi responding point on Fig 3 4 


P ! - a -a 

V 4 b 

r 

£ 

(6) 

6 

, p 

' 2 (a + y) 

0) 

X fa IT 
n V F ’ V a 


(2) 

7 

n 

(i) 

p\V-aW 2 - 

2yW 2 

Kw 

(2) 

8 

— ,X T 

T, ’ 

(4) 

PW - aW 2 - 

~,w 

(2) 

Q 

i pr, 

(4) 

2 yW 2 


T , ‘ 2 { aT , +r) 

x ,W 
nW 


(2 or 3) 




1 , w 


(1 or 2) 





Now if we relax the condition that there is only one vehicle,, then the m opi and V° p! will occur at 
one of the fust seven points only 


3.3.3 Optimisation with respect to n 


Foi the sake of simplicity we will consider the problem when the condition of one vehicle for 
each retailei is lelaxed i e the constraint given by Eqn 3 17 is not effective 
To examine the effect of change of value of n on the optimal solution, we consider the 
following thiee cases 


Case 1 /»* >: 1 




Case II m* > 
Case III 

Case 1 ) When m * > 


> W 


a 


m*< 1 


- < W 


V a 

For small values of //, m*V* < X/n and hence this will remain the optimal to the constrained 
pi oblem I his will continue as n is less than n } ~ 

At this value of n m*V*^X/ni and foi n largei than this value /»*K* > X/n, on the 
boundary and (m* V*) will be outside the feasible region 




Figure 3,5. Optimal m and V when unconstrained optimal is the optimal 
Beyond this point, we set mV-X/fi in Eqn 3 21 to get 


(3 30) 

(3 31) 
a 

As m* and V* is outside the feasible region, the optimal solution is at one of the boundary i e 
point 4 oi point 3 Fuithei as theie is only qne maximum point in the region m £ 0,V> 0, the 
value of optimal solution on the surface along V ()pl = (3/(2(a+ym)) will decrease in both the 
dueclion foi the point 1 i e Z, } < Z 2 £ Z 3 









Further Z 3 > Z2 ( as Z 3 is optimal on the boundaty mnV=X) and hence Z 3 > Z4 or Z 3 is the 
optimal point on the boundaiy 

'I he present values of///* and K* aie maintained till it reaches the boundary m - 1 There the 
optimal shifts to {], } This slnft takes place at the second breakpoint which is achieved at 



Figuie 3,6 : Optimal 111 and V when Ovei supply Constraint is Effective 
1 hus the vauation of F u/ „ and m opl with n for a given retailer is as given below fig 3 7 



Figure 3.7: Variation of m 0 ,,t and V op t with n 


i 



Case 11) Now we considet the case wlien W < . — I his case involves the vehicle capacity 

l a 


also 

When n < 


2 y\VX 


then the oversupply constraint is always satisfied So we further 


[pw - aiv 2 - ;•) 

subdivide this case about this point 

, , 17l ^ 2 yWX 

i) When /?< t r 

[pW-aW 1 - /•) 

If we substitute V W in Eqn 3 23 we get the optimal value for this condition as 
^pW-aW 1 -b ' 




2yW 4 


(3 32) 


This value of nh* always lies befoie the inter section of curves V=W and — — — - as 

2 [a + ym) 



The value of the objective function is optimal at point l in figure 3 6 without the constraint, so 

B 

when the constraint is introduced we will tiv to move along the curve V-~ y r As the 

2(a + yrn) 

function is decreasing as we move further fiom point 1 we try to be as close as possible to 
point 1 on this cuive i e at point 2, but as point 3 is always better than point 2 so 3 is the 
optimal point 

so this value will give the optimal value 

u) When n > 7— — r then * (3.33) 

(pW~ccW 2 -F) n\V 

This solution icmains optimal until the suiface the reaches the boundary as defined by m = 1 
This is achieved at n ~ X/W Beyond this point 1 and V*=X/n 
Thus the variation of V opt and m opl with n is as given below 





Figure 3.8 : Vacation of m op i and V opt with n 

Case HI 

Next we considet the case /»*"' ( 

Again we fuithei subdivide this case also according to the vehicle capacity The vehicle 
capacity constiaint becomes effective only when it cuts the line in =1 below the point where m 

8 . 8 

~ 1 intei sects with V~~, r i e r 

2 [a + ym) 2[a + y) 

a 

When W is greatei than this value then w* = 1 and — r but this value should also 

2{a + y) 


2X(a + y) . _ 

satisfy the ovei supply constiaint so if n > — then r*- X/n 

3W - aW 2 - F 

When W is lessei than this value then V %st W and m 2 *= 

iyw 

P 

As explained earliet this value of lies before the intersection of the curve V~ ~ “j" 


and V~ W but this w* should also satisfy m £ 1 constraint hence if m 2 *^ 1 then ni op( -m 2 * and 
the bieakpoint is X/Wm op , This value of («* V*) remains till the surface nmV-X reaches m - 
1 This is achieved at n *X/W Beyond this point m*= i and V*=X/n When m 2 * £ 1 then the 
fust bieakpoint is missing while the second breakpoint lemains the same 



3.3 4 Solution Method for Problem P 

In the section 3 3 2 we have concluded that for a given value of n for a retailer the optimal 
point will be one of the seven points combination We will refer to these seven points as type 1 
to type 7 solutions 1 he bicakpoints aie being referied to as B1 for the first breakpoint and B2 
as the second one 

We have luitliei divided n into three ranges (I, Bl), (Bl, B2) and (B2, oo) and have identified 
foi each such innge ioi a paiticulai letailei the type of optimal solution (as a function of ??) in 
that lange 


Objective Function 

If m a given lange of ?/, foi a ietailer the optimal solution is either of type 2, 3, 5 or 6 then the 
objective function is a lineai function of n This implies that when all the retailers are such that 
they belong to this set then foi maximising Z we will choose the highest possible value of??. 
When all the ielaileis’ leplenishment is of type two or four or seven then the contribution is of 

B 

the natuic Z~ A~ nC , where A and B aie constants, and when the replenishment is of 

/? 

mixed nature t e some aie being supplied at level independent of n and some at variable level 


then also the naluie of conti ibution is again of type Z = A ~nC , of course with different 

n 

values of constants 


Thus we have again 2 fa 1 pairs of bieakpoints as in the previous case Almost the same 
proccdute is again employed to find the optimal decision variables These breakpoints are 
at ranged in ascending mannei and then paiis are made out of the adjacent numbers Now 
between these tangos we have optimal value of /?? and V and we have to find the optimal ?? 
These tanges will confoim to one of the five cases discussed below 


Case 1, 

When all the values of m/s and V/s aie both constant as m type (1,3, 5, 6) then, as explained 
cailiei, to maximise the objective function value we will try to increase n as large as possible 
i e the uppei limit of the lange will be chosen as the optimal value of?? 



Case 2 


VVlicn the optimal s and V } 's aie all following the relationship 


* X j a j 

n V F. 


and V/- — for a g lven range, then on substituting 
\ K . 


these values in the onginal objective function we can get it as a function of n 
7/„) (If IV .v -III W IF 1 w, fF hv 


W FJjp±^ F . 

> ' 21 / ' pj P p, 2^,’V a, 


b. 1 


J-7 +E-V, +ht 'L x J 1 -nc- 


J J 1 


i r i f * / v (/ j - /V.) k (p - x ) 

Wiiting the above equation in a simplified foim we get 
Z(n)- A - nC - — 


(3 35) 


(3 36) 
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(3,37) 


- ( < 0 ) hence the fn st dei i vati ve gives the maxima (3 3 8) 


(3 39) 


Case 3 


When all the retailci s have their optimal m and V according to the relationship m = —— and 

n\V 

\ 1} , then substituting these values m the original objective function we get 


m 


<">P -*£(>-#, -ftv, 


Jl * b X * 


; i 


, i w, 


-«C- 


U » 1 J*\ l 


(3 40 ) 


This is almost the same as in above case and by similar analysis we get /;*= where 


IT 


S(1 - #>,*, 

W~' * ;=l P 


Now we check the lange foi which this case is true 

If /;* > uppei limit of the lange then n opt = upper limit 

If //*< lowet limit of the lange then n opt - lower limit else n^n*. 


Case 4 

When all the lelaileis have then optimal m and V according to the relationship /«= 1 and 
V X/n then substituting these values in the oiiginal objective function we get again a function 

B 

of the naluic Z(n)- A -n(' and the n* is smulaily obtained as in the above case 

n 

Case 5 

When the leplenishmcnt is of mixed nature, i e some are at constant level and some are at full 
replenishment level, is being considcied in this case 

Let the conti ibulion fiom the constant replenishments be nD+B*> then following the above 
analysis the objective function can be expiessed as 

Z = A vdn-C)-- 


(3 41 ) 


When JM' >0 then the objective function value increases with increase in ti so we 
choose the uppei limit of // m this case 

When (/)-f)< 0 then again following the above analysis we can easily obtain 


n 


B 


C - D 


(3 42) 


Now we check (he uinge foi which this case is tiue 

If//* > uppei limit of the langc then n op( - upper limit 

If //*< lowei limit of the lange then n Q?t - lower limit else n opf =n* 


3.4 Algorithm 

First, the ictailerwise subproblein is solved independently and then the original problem is 
dealt with 
1 Subproblem 

a) Calculate a J} p, and Yj 


a,= (l-ft,) 




IP 2 P, 


2X 


, f>j=R +s-fj- hi 7), and 


i J 


/ / 


l 2 PX, ) 


hi) l- ft, 


' P- X,' 
XlPX^j 


hp 


b) Calculate K ; * : 


1 


/»,*-• 


'P 2 a 
4 F 


- a 


a 


» "V 


c) Calculate the breakpoints with the help of the decision table 3 1 given below 


II Problem 

a) Combine the an ays of breakpoints and then arrange the series in an 
ascending maimer 

b) Make paiis out of the above aira.iged list from each consecutive numbers to 

constiuct ranges 

c) Check the type of eacli range 

d) Calculate the >/„„ and the conesponding m of „ V opt and objective function 


i/nli ta (Vat* onnlr I 



e) I lie lange giving the best objective function value is chosen as the optimal 
taiige and the decision vatiables for this range are the optimal decision 
variables 

When m„* > l 




























3.5 Case When m is Required to be Integer 

Unconsti amed Case 

In the case in* and V* is within the feasible region Further as the function is 
concave on the sui face with respect to m 
V uj»t ~ (l/(2(atym)) 

1 lencc the optimal solution will be L/n*J or and V op f= p/(2(a+yU/* J) 
ot V 0| ,r (3/(2(a 1-yT m*"|) 

In case [ /??*] will fall outside the feasible region, optimal point will be on the boundary 
Const! amed Case 

As the function is concave along the boundary, the optimal integer solution will be 
rounded up oi down value of/;/* and appropimtely adjusted V op , 

3.6 Complexity Analysis 

Let k indicate the numbet of retaileis In the first part, SP, each of the five step requires a 
constant numbet of operations (addition, subtraction, comparison, division etc ) so the total 
numbet of operations foi one retailer is a constant, say L This procedure is repeated k number 
of limes thus the numbei of opeiations i equired is Lk assuming that the square root operation 
takes a constant amount of time i e opeiation time is independent of problem size 
For the pioblem P we have to fust soit the 2k langes which requires {2k log 2k) steps Now 
for each imige we make k evaluations to find the optimal decision variables This step requires 
2k 2 steps 

Thus the pioblem tequiics Lk \ 2k log 2k \ 2k 2 steps Hence the complexity of the above 
algorithm is 0(k 2 )i e it is a polynomial function of number of letaileis 



3.7 Flow Chart 



— 1 

Read the problem 
data 



. _ _ 

Foi each 
calculate the set 

retailer j 
of breakpoints 

i 


Combine the set of breakpoints and 
then soit in an ascending series 


Artange the seues to 
foun langes 


Calculate the n op i fot each range 
according to the range type 


~ — 

Select the decision variables from 

the lange with highest objective 
function value as the optimal 
decision variables 


Figuic 3.9 : Flow Chait for the Algorithm 


3.8 Simultaneous Production-Multiple Replenishment to each 
Retailer 


Wo consider a situation where p.oduction is taking place throughout the production cycle and 
shipments to a retail* mo made at equally spaced intervals throughout the cycle It may be 









mg losI Now we will foimulatc the net profit 


shipment to minimise the iclaiiei invenUny Hold 
function for this case 

3.8.1 Problem Formulation 

A) Ptoduction related cost 

I he following figuies tepiesents inventory status diagram at the production site The 
relevant costs aie cost of production setup and inventory costs 

1 ) Pi eduction set up cost Thei e are n production setups in a unit time As the cost of 

set up is C\ the total pi oduction set up cost per unit time is nC (3 43) 

2) Inventory holding cost I he total inventory holding cost in one period can be 
computed by mcasuiing the aiea of the above figure If the inventory holding cost is hp then 

the total mventoiy cost 

= f (3 44) 

B) Transportation related cost 

1) Cost of transportation 

A fixed fi eight charge Vj and a variable cost per unit transported is incurred 
depending upon the i etailer Thus the total transportation cost is 

»[£/** 1 * (3 45 ) 

r i j=i 

2) Transportation inventory holding cost 

If the tiansportation time is T } for the f retailer and the transportation 
invenloiy holding cost is hi then the total inventory holding cost during transportation over the 

k 

unit time period is n hi (3 46) 

j-i 


, 1 k nijV 

pei unit penod for this is n hp — V ~ 

2 U P 



0 Utivomio and cost at the letailei situ 


1) Revenue 


1 lie total amount of demand satisfied in one period is n V } The total 

/=> 

k 

levenuc genciatcd is n R (3 47) 

2) Inventory holding cost 

Hie following figuies repiesents the stock position at any time at the retailer 

outlet 



Figure 3.10 ; Inventory Status Diagram at Production Site 



Figure 3.1 1 Inventoiy Status Diagtam at Retailer Site 



If the pci unit holding cost at the letailei's site is ht Jy the total inventory holding cost per unit 
time pei lod is 


n 


1 * 


(3 48) 


3) Shot tage cost 


3 he numbei of lost sales in one period at a retailer j is ( X r n WjVj) If the per 
unit shoitage cost is s then the shortage cost per unit time period is 

s jr(Xj-iwVj) (3 49) 

r \ 


The piofit function is obtained by subtracting all costs from the revenue obtained The profit is 
a function of//, ///, and J^’s Let us tepresent this function by Z\ Then this problem OP can be 
stated as following 


Pi oblciii OP 


z( ii, nij, v,)= ii r Y J ,n ’ v ) '»/? *■ T, n, if> y j }~ nht Yj mV,T> 

r i 1 1 j= i ;=* 

-S jr{X-imV,) ( 35 °) 


1 J, 

x, 


j*' 


(3 51) 


We will maximise Z to obtain the optimal policy The feasible policy has to further satisfy the 

following constants 

1) Supply is less than oi equal to demand 

n nijVj *Xj 

2) 3 ho volume shipped is less than or equal to the vehicle capacity 

VjlW s ^ 3 52) 

3) Only one vehicle is allotted per station so the production cycle must be so spaced 

that the vehicle can telurn to take the new replenishment 

(3 53) 


n mint'”] 
h 

4) non negativity consttaint 
/?, ntj Vj £ 0 


.(3 54) 


3.0.2 Solution Methodology 

The eqn 3 50 can be wutten as a function of n as 




t{R - f, - hit ; h - 1 - t"j F j 


c 


n 


k 

'L{Pj»' i v , 


-<2X 


J = 1 

(3 55) 




whci e fij -R-fj -hi 1) \ s and ot, = — + 

2P 2X } 

subject to the constraints as given by Eqn 3 51 to 3 54 


It can be obseived that foi a selected value of n the objective function can be written as 

A 

Z(V)”^T J ) huithei theie aie no constiaint binding the different decision variables together 

j i 

so the ptoblem can be decomposed in k independent subproblems, each of the given type, 


Subproblem OSP 

The i educed /'subpioblem is 

max Z^m,! 7 ) * m{fiV-aV 2 -F) , (3 56) 

l ilts function indicates that foi the const! aints given above, the function will try to be on the 
suiface nmV=X Substituting this value of m in eqn 3 56 we obtain 


Zc~\Xp-aV 

n\ 



(%. FX 

L _ _ a -| __ 

(V V 2 

d% 2 FX 
3V 1 ~ K 3 



( <0 hence the function is concave) 


and m*= 



(3 57) 
(3 58) 
(3 59) 
(3 60) 


awn 


,n3 .RARK 



1 Ins value o( ///* and V* will lemain optimal till tins curve ( nmV = X) touches the boundary at 
//j — 1 1 liis shill; takes place at the value u — XIV* from where tit *— 1 and V * ~ X/n 
If the value of V* obtained in Eqn 3 60 is gieater than W then V*~W and m*=X/nW and 
sinnlaily as the above case this value would remain till it reaches m = 1 boundary This 
intei section takes place at tr X/W fiom wheie///*=l and V*=X/n 

Now the value obtained (oi m * in Eqn 3 60 can be below 1 which is not permissible For this 
case we would try to lemain on the curve wnV—X and reach upto in — 1 on this curve So the 
optimal will be 1 and V* - X/n and if this value is also greater than W then we set V*—W 

Problem OP 

l hus fiom the above analysis we can get the in* and V * at any value of n for each retailer 
I he set of possible m* and V * aie 


1 

* i^ a 

ti>™ J — ~ 

V F 

1 

ii 

» £ 

2 

w*- X/nW 

j y*=w 

3 

nt*~\ 

V*-X/n 

4 

m*~\ 

v*~ w 


This value shifts to othei kind at a bieakpoint The bieakpoints have also been discussed 
above Fit st the langes aie detei mined and then the breakpoints are arranged in an ascending 
mantlet to obtain continuos function between the breakpoints The solution methodology then 
follows the eailiei two cases Fiist the range type is determined and then the local optimal is 
calculated for that tange Fiom the peimissible set of //;* and K* it is clear that the overall 
objective function will be of the kind A-nOB/n The rest of solution methodology is exactly 


the same 


3.9 Example 


I o lllusti ate the above analysis we consider a numerical example with the following data 
C~ $ 1 000 000 pei pioduction set up, R = $200 per unit, P = 30 000 units per year, s = $600 
pci unit, hp = $100 pci unit pet year, ht — $100 per unit per year 
Data assumed to be common foi all letaileis 

(hi = $200 pei unit per yeai, W = 2000 units, f= $2 per unit, F = $100 000 per transportation) 
Xi= 12 000 units T\ - 0 08 yeai 

X 2 = 6 000 units T 2 = 0 05 year 

Xt -= 9 000 units 1 3 = 0 10 year 

Solution 

Initially the i clailct wise subpioblem is solved and the following values are obtained 



Retailer 1 

Retailer 2 

Retailer 3 

a 

0 0125 

0 0125 

0 0167 

IS 

o 

<3\ 

t" 

793 

788 

Y 

0 002 

0 0067 

0 00369 

131 

0 0839 

0 115 

0 099 

B2 

6 00 

3 00 

4 50 

m* 

63 57 

19 05 

37 0 

v* 

2828 

L — — 

2828 

2447 


Table 3.2 Subproblem Calculation 


The breakpoints aic then arranged iit an ascending manner to get piecewtse continuous 
function Between these breakpoints, the optimal value of n is determined and the 
corresponding m„„, and V 0| „ are calculated with the help of the decision table 

1 1re value of lire objective function is calculated for each range A summary of the results for 
the integiatcd pioblcm is given below in fable 3 3 




















Range 

Optimal n 

Objective function 

0-00839 

0 0839 

-7806154 

0 0839 -0 099 

0 099 

-6392449 

0 099-0 115 

0 115 

-5161844 

0 U5 -3 00 

0 956 

456276 8 

3 00 - 4 50 

3 00 

-499999 

4 50- 6 00 

.... — - 

4 50 

-3983332 


Table 3.3: Problem Calculation 


As cieaily evident Horn the above table, the highest objective 
The couespondmg decision vaiiables aie 
mi "6 V j = 2000 

ma - 3 V? “ 2000 

tm-5 " 1881 


function is obtained at n-0 956 




















Chapter 4 
CONCLUSION 
4.1 Conclusion 

I he present wotk dealt with the pioblem of obtaining optimum transportation policy for a 
pioduction distubution system with one production centre and multiple retailers by analysing 
the trade-offs between production, inventory and transportation costs 

Initially the situation when only one replenishment is being made to a retailer in a single 
pioduction cycle was consideied I he pioblem was solved by decomposing it to retatlerwise 
subpiobiem to find the optimal leplemshment volume for a retailer for a given number of 
pioduction set-ups The results indicated that the objective function was a piecewise 
continuous function of numbci of set-ups with the discontinuities at the points where 
replenishment policy changes foi a retailei Between all the bieakpoints the optimal number of 
pioduction set-ups and the couesponding objective function value is calculated and the range 
which gives the best objective function value gives the optimal decision variables 

Latei the case with multiple replenishment was analysed For this situation two types of 
policies single set for each i elailer and multiple set-up for each retailer were analysed In both 
the cases, the pioblem was fust formulated and then solved on the lines of the single 
replenishment case Fiist it was decomposed into i etallei wise subproblem and then solved 
independently Finally the subptoblems weie integrated to get the overall objective function 
which was then solved by dividing it into continues ranges as in the previous case 

The algorithms ptescnled to solve the pioblems have the complexity of 0(^ ), where k is the 
number of t he iclailers i c, it is a polynomial function oi the number of retailers under the 
assumption that square loot opetation of a number can be determined in constant time 


4.2 Scope for future work 

The piesent work can be developed in the following areas 

1 r I he case when only when pioduct is being supplied with only one production centre has 
been consideied This model can be generalised by formulating the problem for multiple 
oi ignis with multiple pioducts 

2 The demand was assumed to be deteimmistic Stochastic demands may be an interesting 
ateafoi fuither leseaich 

3 I he pi oblem can also be analysed foi the condition where the production rate is lower than 
the total demand rate 

4 Integi ation of vehicle routing component will also pose a challenging work 


Appendix 


A.t. 

Aiea of the diagiam lepresenting the inventory status at production site 



< > 


t = V/X 1 

lnventoiy Status Dmgiam at Production Site for Case II 


I he aiea can be divided into three parts the triangles, the rectangles starting from the 
point where the pioduction has stopped and the third part is of the rectangles that fail 
below the tiiangles 
Aiea of fust part 


Nuinbci of triangles = 


mX 

~¥ 


Let the number be denoted by b 


Aiea of tiiangles- 


1 PV 2 b 

2 X 2 


Ai ea of the second pai t 


Aiea of a slab ~ 


PV 
{ X 


V 


v_ 

X 


Numbci of slabs - b(b-l)/2 

r/2 

Hence the area = •— j(P ~ X)b(b - l) 
2X 


Aiea of t h e third pai t 


The numbei of blocks aftei the production has stopped is 


m(in + 1) 

_ 


-b 


t Ience the ai ea of this pai t is - 




1 he total area is calculated by adding all the three parts obtained above 
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